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SUMMARY 


This  report  addresses  the  question  of  the  coupling  between  the  ionosphere 
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containing  hot  plasma. 

Previous  analyses  of  the  S3-3  satellite  data  indicate  that  the  magnetic- 
field-aligned  current  sheet  above  the  aurorae  generally  becomes  unstable, 
generating  electrostatic  wave  turbulence  and  high  potential  differences  along 
the  magnetic  field.  Here,  we  present  information  on  the  conditions  in  the 
unstable  current  region  inferred  from  the  analysis  of  data  obtained  along  an 
unusual  trajectory  of  the  satellite  wherein  the  satellite  remained  in  and 
near  the  current  sheet  for  an  interval  of  about  six  hours  in  magnetic  local 
time.  Potential  differences  along  the  magnetic  field,  above  and  below  the 
satellite  (altitude  -7800  km)  are  inferred  from  the  energy  and  pitch-angle 
distributions  of  the  electrons.  These  potential  differences,  which  exceeded 
10  kV  at  times,  were  generally  higher  during  the  local  evening  than  during 
the  local  afternoon.  They  were  also  quite  variable  along  the  satellite 
trajectory,  indicating  latitudinal  fluctuations  of  the  current  sheet  with 
periods  of  2-7  minutes. 

Upward-flowing  H+  and  0+  ions  appeared  whenever  the  potential  differences 
exceeded  about  .5  kV,  which  corresponds  to  the  lowest  energy  threshold  of  the 
detectors.  However,  the  pitch-angle  distributions  of  the  ions  imply  that  the 
ions  were  accelerated  not  only  along  the  magnetic  field  by  the  potential 
differences,  but  also  perpendicularly  to  the  field.  The  transverse  accelera¬ 
tion  of  the  ions  most  probably  occurred  through  interactions  with  electro¬ 
static  ion  cyclotron  (EIC)  waves,  which  are  an  expected  product  of  the  current- 
driven  instability,  and  which  have  been  observed  in  the  current-sheet  region. 

A  Monte  Carlo  program  was  used  to  compute  energy  and  pitch-angle  distri¬ 
butions  of  ions,  taking  into  account  the  effects  of  the  waves,  the  potential 
difference  along  the  magnetic  field,  and  the  divergence  of  the  magnetic  field. 
Two  measurements  of  the  ion  distribution  were  reconstructed  by  choosing  a 
transverse-heating  rate  appropriate  for  measurements  of  the  wave  power  spectrum, 
a  potential  difference  below  the  satellite  equal  to  that  inferred  from  the 
electron  distribution,  a  source  base  2000  km  to  4000  km  below  the  satellite, 
and  a  source  intensity  that  decreases  exponentially  at  altitudes  above  the 
base  with  a  scale  height  of  600  km  to  1000  km.  An  ion  distribution,  measured 
when  the  potential  difference  above  the  satellite  was  high,  was  reconstructed 
with  parameters  similar  to  the  above,  except  with  a  source  intensity  that 
increased  with  altitude  above  the  base.  The  parameters  used  to  fit  the 
measured  distributions  are  consistent  with  anomalous-resistivity  theory. 


1 


Recommendations  are  made  for  further  work  that  is  necessary  in  order 
to  assess  the  effects  of  the  current-driven, instability  on  the  environmental 
conditions  resulting  from  a  nuclear  explosion. 
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SECTION  I 


INTRODUCTION 


The  electrical  coupling  of  the  ionosphere  and  magnetosphere  strongly 
affects  several  components  of  the  environment  resulting  from  a  high-altitude 
nuclear  explosion.  The  most  important  of  these  are  1)  the  plasma  irregularities, 
which  degrade  communications  and  radar  systems,  2)  the  emissions  due  to  energy 
and  debris  deposition  in  the  atmosphere,  which  enhance  the  backgrounds  of 
optical/IR  systems,  3)  the  distribution  of  the  energetic  ionized  debris  and 
air  in  the  magnetosphere,  which  erode  satellite  surfaces,  and  4)  the  distri¬ 
bution  of  the  trapped  relativistic  electrons,  which  degrade  satellite 
components.  In  all  the  work  that  has  been  done  to  date  on  the  delineation  of 
these  environmental  components,  the  classical  values  of  the  conductivity  along 
the  magnetic  field,  based  on  binary  collisions,  have  been  used.  Above  the  F- 
layer  of  the  ionosphere  such  conductivities  are  so  low  that  the  magnetic  field 
lines  have  been  regarded  to  be  equipotent ials .  All  this  work  now  needs  to  be 
reassessed  in  view  of  the  S3-3  Satellite  data,  which  unequivocally  imply 
the  presence  of  high  electrostatic  potential  differences  along  the  magnetic 
field . 


The  data  reveal  that  these  potential  differences  extend  to  tens  of 
kilovolts  (Refs.  1,  2,  3)  and  are  nearly  always  present  in  a  narrow  magnetic 
shell,  of  thickness  equal  to  a  few  degrees  in  latitude,  connected  to  the 
auroral  oval  (Ref.  4).  Specifically,  the  elevated  potentials  occur  in  the 
region  of  the  upward  flowing,  magnetic-field-aligned  currents  described  by 
lijima  and  Potemra  (Ref.  3).  They  are  generally  below  the  altitude  of  the 
S3-3  satellite  (~1  earth  radius),  but  occasionally  extend  to  well  above  the 
satellite.  The  potential  gradients  are  predominantly  in  a  direction  such 
that  electrons  are  accelerated  downward  and  ions  upward.  Indeed,  ions 
(principally  H+  and  0+)  in  the  keV  range  are  often  observed  moving  upward, 
generally  closely  aligned  with  the  magnetic  field,  when  these  potential 
differences  occur  below  the  satellite  (Ref.  4,  6,  and  7). 


The  pitch-angle  distributions  of  the  ions,  however,  are  much  broader 
than  expected  if  ions  of  the  thermal  population  were  simply  accelerated  along 
the  magnetic  field.  The  ions  must  also  have  been  accelerated  in  the  direction 
perpendicular  to  the  magnetic  field.  Acceleration  by  waves  is  suggested. 
Turbulent  electrostatic  plasma  waves,  including  electrostatic  ion-cyclotron 
(EIC)  waves,  are  nearly  always  observed  in  the  region  where  the  potentials 
occur  (Refs.  8  and  9).  It  is  therefore  likely  that  the  ions  acquired  their 
transverse  energies  by  resonating  with  the  EIC  waves.  In  the  adjacent  regions 
where  the  field-aligned  current  flows  downward  (Ref.  3),  coherent  EIC  waves 
are  observed.  There,  the  pitch-angle  distributions  of  energetic  H+,  He+  end 
0+  ions  are  found  to  be  conical,  i.e.,  peaked  at  pitch-angles  between  0  and 
90°  (Ref.  1°).  These  distributions  indicate  that  the  ions  are  rapidly 
accelerated,  within  a  short  distance  along  the  magnetic  field,  by  the  coherent 
EIC  waves. 


The  satellite  measurements  described  above  indicate  that  a  current- 
driven  instability,  generating  electrostatic  waves,  may  provide  the  mechanism 
necessary  to  support  the  field-aligned  potential  differences  in  the  absence 
of  collisions.  Three  theories  of  such  instabilities  have  been  proposed:  The 
electric  double  layer  (Refs.  11  and  12),  the  blique  electrostatic  shock 
(Ref.  13),  and  anomalous  resistivity  (Ref.  14  and  15).  Some  features  of  each 
theory,  at  different  times,  seem  to  be  implied  by  the  data.  The  Kindel  and 
Kennel  theory  (Ref.  14),  however,  appears  to  explain  the  majority  of  the 
observations.  Their  theory,  for  example,  predicts  the  following: 

(i)  A  field-aligned  current  will  become  unstable  above  the 
F-layer  of  the  ionosphere  in  the  density  range  2-20pA/m^. 

(Such  densities  have  been  confirmed  by  measurements.) 

(ii)  The  instability  will  produce  electrostatic-wave  turbulence 
at  drift  velocities  that  are  high  compared  with  the  electron 
thermal  velocity.  (The  turbulence  observed  in  the  upward- 
flowing-current  region  may  therefore  be  due  to  the  high 
drift  velocity  of  the  downward-moving  current-carrying 
electrons) . 

(iii)  EIC  waves  will  be  destabilized  over  a  wide  range  of  electron 
to  ion  temperatures  (0.2  to  8)  at  the  lowest  current 
densities  and  critical-drift  velocities.  (The  EIC  waves  in 
the  absence  of  turbulence  observed  in  the  downward-f lowing 
current  region  may  therefore  be  due  to  the  low  drift  velocity 
of  the  upward-moving  electrons.) 

(iv)  In  a  plasma  containing  two  ionic  constituents,  the  heavier 
ion  cyclotron  waves  will  be  more  easily  destabilized,  even  if 
the  heavier-ion  abundance  is  as  low  as  about  10%.  (Although 
the  low-frequency  He+  and  0+  cyclotron  waves  are  difficult  to 
identify  because  of  the  noise,  their  effects,  i.e.,  the 
conical  distributions  of  the  He+  and  0+  ions,  have  been 
observed . ) 

In  the  context  of  the  Kindel  and  Kennel  theory,  the  potential  differences 
are  due  to  anomalous  resistivity.  The  wave  turbulence  provides  the  resistive 
drag  that  inhibits  the  directed  motion  of  the  electrons  along  the  magnetic 
field  and  thus  allows  the  potentials  to  build  up.  However,  in  order  to 
appreciably  retard  the  electron  motion,  it  is  necessary  for  the  electric 
fields  of  the  waves  to  be  larger  than  the  field  aligned  electric  field  arising 
from  the  potential  (Ref.  16).  This  requires  that  the  scale  of  the  observed 
potentials  be  large:  The  measured  wave  fields  are  less  than  about  10  mV/m 
(Ref.  9)  and  the  potential  differences  are  generally  larger  than  1  kV,  hence, 
the  scale  of  the  potential,  and  the  extent  of  the  wave  turbulence,  must  be 
of  the  order  of  10^V/ 10~3v/m  =  10^m,  or  1000  km,  along  the  magnetic  field. 

Quite  often,  however,  electric  fields  as  high  as  about  0.5V/m  are  measured 
(Ref.  8),  implying  shorter  scales  of  the  potentials.  Hence,  double  layers  or 
electrostatic  shocks  might  also  occur  in  the  field-aligned  current  regions. 


It  is  clear  that  much  more  information  on  these  natural  processes 
must  be  obtained  before  we  can  hope  to  predict  the  nuclear-induced  environ¬ 
ment.  In  particular,  information  is  needed  on  the  distribution  (location 
and  extent)  of  the  potential  along  the  magnetic  field  and  on  the  nature  and 
distribution  of  the  wave  turbulence. 

We  have  been  studying  the  potential  distributions  by  analyzing  the 
pitch-angle  and  energy  distributions  of  the  electrons  and  ions  obtained 
with  the  Lockheed  experiment  (Ref.  6)  on  the  S3-3  satellite.  In  our  previous 
DNA  report  (Ref.  1)  we  described  the  "inverted-V"  potential  structure  that 
was  inferred  from  the  observation  while  the  satellite  passed  almost  normally 
through  the  current  sheet.  In  this  report  we  present  further  information, 
based  on  the  satellite  measurements,  on  the  properties  of  the  instability  and 
the  conditions  under  which  it  occurs. 

In  Section  II  we  point  out  the  similarity  of  physical  processes 
associated  with  the  convection  of  field  lines  over  the  polar  cap  to  the 
processes  resulting  from  the  expansion  of  the  hot  plasma  produced  by  a 
nuclear  explosion  at  high  altitudes.  This  analogy  clearly  reveals  that 
the  auroral-zone  observations  discussed  above  are  directly  relevant  to  the 
ionosphere-magnetosphere  coupling  resulting  from  a  nuclear  burst.  In 
Section  III  we  present  the  Lockheed  data  obtained  along  an  unusual  trajectory 
wherein  the  satellite  remained  in  and  near  the  current  sheet,  within  a  few 
degrees  of  invariant  latitude,  over  a  local  time  interval  of  about  6  hours; 
and  we  discuss  the  conditions  in  that  region  inferred  from  an  analysis  of 
the  data.  In  Section  IV  we  discuss  the  theory  of  the  pitch-angle  diffusion 
of  ions  due  to  interactions  with  EIC  waves,  in  the  presence  of  an  electric 
field  along  the  magnetic  field,  and  describe  a  Monte  Carlo  calculation  that 
was  used  to  compute  the  ion  distributions.  The  properties  of  the  EIC  wave 
turbulence  inferred  from  the  comparison  of  the  computed  ion  distributions 
with  the  measured  distributions  are  discussed.  The  conclusions  and  recommen¬ 
dations  for  further  work  are  discussed  in  Sections  V  and  VI,  respectively. 


SECTION  II 


RELEVANCE  OF  POLAR-CAP  PROCESSES 
TO  NUCLEAR-BURST-INDUCED  PROCESSES 


1.  CONVECTION  OF  MAGNETIC  TUBE  CONTAINING  HOT  PLASMA 

Following  a  nuclear  explosion  at  high  altitudes,  hot  plasma  expands 
into  a  local  magnetic-tube  of  force,  as  depicted  in  Figure  1.  Initially, 
this  plasma  consists  of  the  bomb  debris,  including  fission  fragments  and 
relativistic  fission  beta-decay  electrons,  and  atmospheric  constituents  which 
have  been  energized  and  ionized  near  the  burst  point.  Owing  to  the  plasma 
pressure  and  the  magnetic-field  configuration,  an  outward  force  acts  on  the 
tube  and  polarizes  the  plasma  in  the  ai imuthal  direction.  As  shown  in  Figure 

1,  a  net  negative  charge  appears  at  the  eastern  surface  of  the  tube,  and  a 

net  positive  charge  appears  at  the  western  surface.  An  electric  field  directed 
toward  the  east,  therefore,  develops  in  the  tube,  and  the  plasma  drifts  out¬ 
ward  with  the  E  X  B  drift  velocity.  In  this  manner  the  debris  and  the  trapped 
relativistic  electrons  become  distributed  to  regions  of  the  magnetosphere  above 
the  magnetic  shell  in  the  vicinity  of  the  burst. 

At  somewhat  later  times,  the  heated  plasma  resulting  from  energy 
deposition  in  the  atmosphere  expands  into  a  local  magnetic  tube  and  polarizes 
either  in  the  direction  discussed  above  or  in  the  opposite  direction, 
depending  on  the  magnitude  of  its  kinetic-energy  density.  In  either  case,  the 
combination  of  the  resulting  plasma  E  X  B  drift  velocity  and  density  gradient 
cause  the  plasma  to  striate.  The  striated,  high-density  plasma  thereby  becomes 
distributed  over  large  distances  in  the  magnetosphere  and  interferes  severely 
with  radio-wave  propagation. 

Both  the  distribution  of  the  debris  at  early  times  and  the  plasma 
striations  at  later  times,  depend  sensitively  on  the  coupling  of  the  ionosphere 
and  magnetosphere.  The  charge  separation  across  the  magnetic  tube  tends  to 
be  neutralized  by  currents  in  the  thermal  plasma  which,  as  shown  in  Figure  1, 
flow  along  magnetic  field  lines  at  high  altitudes  and  across  magnetic  field 
lines  in  the  ionosphere,  where  the  collision  frequency  is  sufficiently  high. 

In  the  applications  mentioned  above,  the  neutralizing  current  has  been  accounted 
for  by  using  the  classical  value  of  the  conductivity.  If  a  potential  difference 
should  develop  in  the  magnetic-field-aligned  current  regions,  due  to  anomalous 
resistivity  or  any  other  cause,  it  would  greatly  alter  the  distributions  of  the 
plasmas . 

2.  CONVECTION  OF  POLAR  CAP  FIELD  LINES 

The  magnetic  field  lines  over  the  polar  cap  are  convected  in  the  antisolar 
direction  by  the  pressure  of  the  solar  wind.  The  motion  induces  a  dawn-to-dusk 
electric  field  over  the  polar  cap,  as  depicted  in  Figure  2.  The  net  positive 
and  negative  charges  are  on  the  "stationary"  field  lines  over  the  aurora  oval, 
the  positive  charge  being  located  on  the  half  of  the  auroral  oval  on  the  dawn 
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ire  1.  Illustration  of  electric  field  across  hot-plasma  tube 
current  tending  to  neutralize  the  excess  charges  at  the  east 
west  boundaries  of  the  tube. 


POLAR-CAP  FIELD  LINE  MOTION 


Figure  2.  Illustration  of  electric  field  over 
polar  cap  and  current  tending  to  neutralize  the 
excess  charges  at  the  dawn  and  dusk  boundaries 
of  the  antisolar  convection  region. 
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side  and  the  negative  charge  on  the  dusk  half  of  the  oval.  The  primary 
currents  that  tend  to  neutralize  the  charges  are  also  shown  in  the  figure. 
They  flow  down  the  field  lines  containing  the  net  positive  charge  on  the 
dawn  side,  across  the  polar  cap  in  the  ionosphere,  and  up  the  field  lines 
containing  the  net  negative  charge.  These  are  the  region  1  currents  identi¬ 
fied  by  lijima  and  Potemra  (Ref.  5).  A  less  intense  region  2  current  system 
is  also  observed  at  high  values  of  Kp.  These  currents  flow  along  field  lines 
adjacent  to  the  region  1  currents  but  at  lower  latitudes,  and  are  directed 
opposite  to  the  region  1  currents. 

The  potential  across  the  polar  cap  is  generally  about  40  kV  and  the 
ionospheric  resistance  over  the  polar  cap  is  about  0.2  ohms.  Hence,  the 
field-aligned  region  1  current  is  about  4  x  loVo.2  =  2x10^  Amperes.  Half 
the  cross  sectional  area  of  the  current  shell  at  the  altitude  of  the  Triad 
satellite  is  approximately, 

A  =  Tr(R£+hq)d  Sin  0  (1) 

where  R^  is  the  earth's  radius,  hs  is  the  altitude  of  the  satellite,  d  is 
the  width  of  the  shell  at  hs ,  and  0  is  the  colatitude  of  the  auroral  zone. 
Taking  hs  =  800  km,  d  =  10  km,  and  0  =  25°,  we  find  that  A  »  9.5xl0^km^, 

The  field-aligned  current  density  is  therefore  2x10^/9.5x10^  =  2,lpA/m  . 

This  value  is  consistent  with  the  current  densities  measured  with  the  Triad 
satellite  (Ref.  5).  It  is  also  within  the  range  of  the  critical  current 
density  given  by  Kindel  and  Kennel  (Ref.  14)  for  the  onset  of  the  current- 
driven  instability.  Moreover,  a  current  density  of  about  lpA/m^  was 
measured  on  the  S3-3  satellite  when  wave  turbulence  and  other  products  of 
the  instability  were  observed  (Ref.  8).  It  is,  therefore,  not  surprising 
that  potential  differences  along  magnetic  field  lines,  as  well  as  magnetic- 
field-aligned  currents,  are  regular  features  of  the  auroral  oval. 

The  field  aligned  currents  along  the  eastern  and  western  boundaries 
of  the  magnetic  tube  containing  the  early-time  plasma,  for  a  Starfish-  or 
USSR1-  type  burst,  is  about  1(PA  (Ref.  17).  Since  the  cross  section  of 
half  the  boundary  just  above  the  E-layer  of  the  ionosphere  is  about  3xl0^m^, 
the  field-aligned  current  density  is  about  30jjA/m^.  The  current  density 
is  the  most  important  parameter  that  determines  whether  the  current  becomes 
unstable  (Ref.  L4) .  Hence,  since  this  value  exceeds  the  critical  value  it 
is  highly  probable  that  the  instability  will  occur  in  the  hot-plasma  tube 
as  it  does  in  the  polar  cap.  Furthermore,  since  the  processes  responsible 
for  the  field-aligned  currents  are  similar,  the  information  that  is  now 
being  obtained  with  the  S3-3  satellite  on  the  nature  and  effects  of  the 
instability  in  the  auroral  zone  should  be  directly  applicable  to  the  plasma 
tube  situation. 
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SECTION  III 

CONDITIONS  IN  UN STABLE- CURRENT  REGION 


1.  S3- 3  SATELLITE  MEASUREMENTS 

The  Lockheed  experiment  (Ref.  6)  on  the  S3-3  satellite  measures  the 
energy  and  pitch-angle  distributions  of  ions  in  the  energy-per-charge  range 
.50  -  J6  keV  and  electrons  in  the  energy  range  70  eV  -  24  keV.  The  measure¬ 
ments  are  obtained  with  three  ion  mass  spectrometers  and  four  magnetic  elec¬ 
tron  spectrometers  mounted  such  that  their  view  directions  are  perpendicular 
to  the  spin  axis  of  the  satellite.  The  spin  rate  of  the  satellite  is  -3.5  RPM 
about  an  axis  perpendicular  to  the  orbital  plane;  hence  a  nearly  complete  fitch- 
angle  scan  is  obtained  in  about  8.7  sec.  The  ion  spectrometers  sample  the  mass- 
per-charge  (M/Q)  distribution  in  the  range  1-30  once  per  second.  Each  ion 
spectrometer  has  4  energy-per-charge  settings  which  are  stepped  every  16  sec. 

The  energy  settings  of  each  ion  spectrometer,  together  with  the  energy  response 
and  geometric  factors  of  the  electron  spectrometers,  are  listed  in  Table  1. 

On  17  January  1977,  the  satellite  remained  in  and  about  the  unstable- 
current  region  for  a  local-time  interval  of  about  6  hours.  The  trajectory  of 
the  satellite  is  shown  in  Figure  3.  There,  the  altitude  and  invar ient  latitude 
are  plotted  against  satellite  spin  number  and  universal  time  in  seconds  (listed 
on  abscissa  at  bottom  of  graph)  and  magnetic  local  time  (listed  at  top  of  graph) . 
The  satellite  spin  number  is  a  useful  parameter  in  the  discussion  of  the  data. 
Note  that  the  satellited  altitude  and  invariant  latitude  change  very  little 
as  the  satellite  moves  from  21.5  h  to  15.7  h  MLT,  during  a  time  interval  of 
about  0.6  h.  All  the  magnetic  indices  were  low  for  several  days  prior  to  this 
satellite  pass.  At  the  time  of  the  measurements  the  three-hourly  Kp  index 
was  2. 


The  data  obtained  on  this  pass  are  shown  in  Figure  4.  At  the  bottom 
of  the  chart  are  shown  the  universal  time  (SYST) ,  longitude,  latitude,  alti¬ 
tude  in  km,  invariant  latitude  (ILA) ,  and  magnetic  local  time  in  h.  The  four 
lowest  panels  show  the  logarithm  of  the  counts  per  half-second  of  the  electron 
spectrometers.  The  panel  labeled  PITCH  shows  the  pitch  angle  of  the  measured 
particles.  Here,  the  pitch  angle  of  0  denotes  particles  moving  directly  down 
the  field  lines.  The  next  four  panels  show  the  logarithm  of  the  sum  of  the 
counts  per  second  of  the  3  mass  spectrometers  for  ions  of  M/Q  =  1,  2,  4,  and 
16  respectively.  A  code  designating  the  energy  steps  of  the  mass  spectro¬ 
meters  is  in  the  top  panel.  The  numbers  on  the  pitch  angle  trace  denote  the 
spin  number  of  the  satellite.  We  will  refer  to  the  ions  which  have  pitch- 
angle  distributions  peaked  at  a=180  ,  as  field-aligned  ions,  although  they 
are  not  strictly  field  aligned  as  discussed  in  the  introduction. 
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Table  1.  Detector  characteristics 


Detector 

Particle 

Energy,  keV 

GDE ,  cm 

sr  keV 

CMEA 

Electrons 

0.07  -  0.24 

1.2  x 

10-6 

CMEB 

Electrons 

0.35  -  1.1 

6.5  x 

icf6 

CMEC 

Electrons 

1.6  -  5.0 

1.9  x 

10-5 

CMED 

Electrons 

7.3  -  24 

6.5  x 

icf5 

Energy  per  unit 

charge  (keV) 

Step  1  2 

3 

4 

CXA  1 

Ions 

0.50  0.68 

0.94 

1.28 

CXA  2 

Ions 

1.76  2.4 

3.3 

4.5 

CXA  3 

Ions 

6.2  8.5 

11.6 

16.0 
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MAGNETIC  LOCAL  TIME  (hr) 


Figure  3.  Altitude  and  invariant  latitude  of  satellite  versus 
satellite  spin  number,  universal  time,  and  magnetic  local  time. 
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While  examining  the  data,  note  the  following  types  of  features: 


(i)  Field-aligned  and  conical  pitch-angle  distributions  of 
H+(M/Q=1),  0+(M/Q)=16,  and  occasionally  h£(M/Q=4)  are  present. 

For  example,  on  spins  2,  3,  and  4,  the  counting  rates  of  H 
and  0+  ions  near  a=180  are  appreciably  above  the  noise  level. 

(ii)  The  troughs  in  the  CME  counting  rates  centered  at  a=180°  are 
wider  and  deeper  when  the  field-aligned  ions  appear.  This 
widening  and  deepening  of  the  loss  cones  of  the  electron 
distributions  is  due  to  an  electrostatic  potential  difference 
along  the  magnetic  field,  below  the  satellite,  that  accelerates 
electrons  downward  and  ions  upward  (Refs.  1  and  2).  This  effect 
is  more  pronounced  for  the  lower  energy  electrons.  Note,  for 
example,  that  the  electron  loss  cones  are  wider  and  deeper  on 
spin  2  than  they  are  on  spin  1. 

(iii)  The  counting  rates  of  the  lower-energy  detectors  are  sometimes 
symmetrical  about  a=90  ,  while  the  counting  rates  of  the  higher 
energy  electrons  are  high  near  a=0  and  have  local  minima  at  a 

=90  .  Such  distributions  indicate  the  presence  of  electrostatic 
potential  differences  above  the  satellite  that  accelerate  elec¬ 
trons  downward  (Refs.  1  and  2).  The  symmetrical  distributions 
are  due  to  electrons  that  are  locally  trapped  by  the  magnetic 
field  below  the  satellite  and  the  electric  field  above.  Note, 
for  example,  that  on  spins  1,  2,  and  3,  the  counting  rates  of 
the  CMEB  detector  are  symmetrical  about  90°  while  the  CMEC 
counting  rates  have  maxima  near  0°  and  local  minima  at  90°. 

(iv)  The  CMEC  counting  rate  becomes  modulated  when  electron  "spikes"— 
high  fluxes  of  electrons  closely  aligned  with  the  magnetic  field — 
appear.  Such  field-aligned  electrons  are  usually  observed  moving 
simultaneously  upward  and  downward  along  the  magnetic  field,  with 
comparable  fluxes  in  both  directions  (Ref.  18  and  19).  On  spins 
7,  8,  and  9,  the  spikes  are  directed  only  upward,  as  indicated  by 
the  CMEA  and  CMEB  counting  rates.  The  largest  spikes  (both  up 
and  down)  appear  on  spins  47-54. 

2.  ANALYSIS  OF  DATA 

In  references  1  and  2  several  methods  were  used  to  determine  the  poten¬ 
tial  differences  above  and  below  the  satellite  in  order  to  evaluate  the 
effects  of  some  of  the  uncertain  parameters,  such  as  the  shape  of  the  poten¬ 
tial  along  the  magnetic  field,  the  electron  spectrum  in  the  electric-field- 
free  region,  and  interaction  of  the  electrons  with  the  wave  turbulence.  Now 
that  this  groundwork  has  been  done,  the  most  direct  methods  can  be  used  to 
estimate  the  potential  differences.  Accordingly,  in  this  analysis,  the 
potential  differences  below  the  satellite  were  determined  from  the  loss  cone 
widths . 
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It  will  be  recalled  that  at  two  points,  (B  ,$  )  and  (B,<J>),  along  the 
dynamical  trajectory  of  an  electron,  the  energy,  w,  and  pitch  angle,  ot,  of 
the  electron  are  related  by  the  equations  expressing  the  conservation  of  the 
magnetic  moment  of  the  electron  and  the  total  energy,  viz.. 


w 


and 


s  .  2  w  .  2 

—  sm  a  =  —  sm  a 
B  s  B 
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w  =  w  +  e($-$  ) 
s  v  s' 
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(3) 


Here,  B  is  the  magnetic  field  intensity,  $  is  the  electrostatic  potential, 
and  e  is  the  absolute  value  of  the  electron  charge.  The  subscript,  s, 
denotes  the  values  of  the  parameters  at  the  location  of  the  satellite.  Hence, 
if  B^  designates  the  limiting  field  intensity  at  the  "top"  of  the  atmosphere 
where  the  electron  mirrors  before  suffering  collisional  effects,  and  where 


$  ,  the  edge  of  the  loss  cone, 
m*  &  9 


-  a 


is  given  by  the  equation 
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The  potential  towarc  lower  altitudes,  is  assumed  to  increase  monotonically 

from  0  to  its  maximum  value,  $  ,  and  then  to  remain  constant. 

m 

This  equation  was  used  to  determine  $  -  0  .  Generally,  the  CMEB 

detector,  which  has  a  high  counting  rate  and  is  sensitive  to  a  relatively 
narrow  range  of  electron  energies,  was  used  to  determine  a  .  For  this 
detector  the  energy  0.9  keV  was  used  for  w  .  However,  the  value  of  a 
obtained  from  the  CMEC  counting  rate  was  often  used  to  verify  results? 

For  this  detector,  the  energy  1.8  keV  was  used  for  wg  in  Eq .  (4).  The  48- 
term  Jensen  and  Cain  magnetic  field  model  for  1972  (Ref.  20)  was  used  to 
determine  the  values  of  Bg  and  B^  on  the  field  lines  traversed  by  the  satel¬ 
lite  trajectory.  B  was  evaluated  at  the  altitude  of  200  km.  The  ratio  B  /B 
along  the  satellite  trajectory,  as  a  function  of  the  satellite  altitude,  i!  t 
shown  in  Figure  5. 

The  potential  at  the  satellite  and  the  value  of  B^  above  the  satellite 
where  <J>  =  0,  were  estimated  from  the  CMEC  counting  rate  by  assuming  the  primary 
electrons  at  B1  to  be  isotropic  in  the  downward  hemisphere  and  to  have  a  Max¬ 
wellian  distribution  with  a  temperature,  w  ,  equal  to  1  keV.  Such  a  distribu¬ 
tion  is  consistent  with  many  samples  of  the  data  obtained  at  times  when  the 
potential  was  zero.  The  directional  flux  of  these  electrons  at  the  satellite, 

B  ,  from  Liouville's  theorem,  is 
s  * 
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Figure  5.  Magnetic-field  ratio  Bs/Bt  as 
of  satellite  altitude. 
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where,  from  the  conservation  of  the  magnetic  moment, 
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for  w  e<ps*  Here,  is  the  pitch  angle  of  the  electrons  at  B]_.  The 
counting  rate  of  the  CMC  detector  due  to  this  flux  is  given  by  the 
integral , 
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G  =  5.4  x  10  cm  sr  is  the  geometric  factor  of  the  detector,  and  is 
the  larger  of  -  1.6  keV  (the  low-energy  limit  of  the  detector  bandpass) 
or 
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which  follows  from  (6)  for  =  90°.  Note  that  the  flux  given  by  (5)  is 
constant  over  the  pass  band  of  the  detector  in  the  pitch-angle  range  0  i 
ots  .1  otL,  where  a ^  follows  from  (8)  for  *  w^.  Hence,  the  counting  rate 
of  the  detector  should  be  constant  in  that  pitch  angle  interval.  The  situa¬ 
tion  is  illustrated  in  Figure  6  which  depicts  the  response  of  the  detector 
as  a  function  of  pitch  angle,  on  the  contours  of  constant  flux  obtained  by 
transforming  the  isotropic  flux  from  B^,  <f>  «  0  to  Bs,  <f>g.  The  relationship 
between  <J>S  and  B^  at  the  limit  of  the  constant  counting-rate  interval,  given 
by  the  equation 


e4,s  =  \  (1  ~  r  sin\} 

s 


(9) 


is  valid  for  an  isotropic  flux  at  B^ ,  regardless  of  the  energy  distribution. 
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Figure  6.  Illustration  of  transf orma t ion  of  flux  which  is 
isotropic  at  B^,  ♦  =  0  to  Bs ,  cf'g,  where  B^/Bs  =  .77  and 
<j>s  =  .8  kV.  The  flux  contours  between  the  numbered  dots  in 
(a)  correspond  to  the  contours  between  the  correspondingly 
numbered  dots  in  (b) .  Note  that  the  counting  rate  of  the 
detector,  which  responds  to  energies  1.6  to  5.0  keV,  is 
constant  between  a  =  0  to  a  =  a^. 


By  putting  (5)  into  (7)  and  integrating  for  a£ 
ratio  C(aL)/C(90  )  is  found  to  be. 


aL  and  ag  -  90°,  the 


C(°tL)  =  (1  +  1.6)  exp  (-1.6)  -  (1  +  5.0)  exp  (-5.0) 
C(90°)  (1  +  exp  (-w.^)  -  (1  +  5.0)  exp  (-5.0) 


where  w..  is  given  by  Eq.  (8)  for  a  =  90°,  i.e., 
1  ® 


W1  ~  1  -  B- /B 
1  s 


The  ratio  C(otj)/C(900)  is  plotted  against  in  Figure  7.  The  straight  lines 
in  Figure  8  give  the  potential  <f>  as  a  function  of  B^/Bg  for  fixed  values  of 
w^.  By  eliminating  4>s  from  Eqs.  (9)  and  (11),  we  obtain 


°l  ■  ^  t1  ■  ^ w  ■  cl 


The  curves  in  Figure  8  show  the  values  of  as  a  function  of  B^/B^  for  fixed 
values  of  w^ .  Hence,  s  and  B^/Bs  for  values  of  a}>s  <.  wl  *  1.6  keV,  are 
determined  from  the  measured  values  of  aL  and  C(aL)/C(90°)  as  follows:  the 
value  of  w-l  corresponding  to  the  counting-rate  ratios  is  obtained  from  the 
graph  of  Figure  7;  B^/Bs  corresponding  to  and  w^  is  obtained  from  Figure  8 
the  4>s  corresponding  to  the  same  values  of  B^/Bs  and  w^  is  read  off  the  graph 
in  Figure  8. 

For  e<(>s  >  wL  the  peak  of  the  CMEC  counting  rate  is  at  a=0,  and  we  do 
not  have  sufficient  information  from  this  detector  alone  to  determine  <J>  and 
Bl/Bs  separately.  However,  we  can  determine  the  lower  and  upper  limits  of 
these  values.  The  counting  rate  ratio,  C(0°) /C (90°) ,  of  the  CMEC  detector 
now  becomes 

c^qOj  (1  +  e^g)  exp  (-e<t>s)  -  (1  +  5.0)  exp  (-5.0) 

C(90°)  (1  +  Wj_)  exp  (-W;l)  -  (1  +  5.0)  exp  (-5.0)  (12) 

After  substituting  Eq .  (11)  for  Wi ,  this  equation  was  solved  numerically  for 
4>s  for  various  values  of  R  =  C(0°)/C(90o)  and  B^/Bs .  The  results  are  shown 
in  Figure  9,  where  <t>s  is  plotted  against  Bl/Bs  for  fixed  values  of  R.  This 


figure  shows  that  for  R  >  1,  4>s  is  between  1.6  -  5.0  kV,  and  Bj_/Bq  is  between  0 
and  the  limiting  value  (B^/B^^  given  by  the  intercept  of  the  appropriate  R 
curve  at  the  abscissa.  If  <J>S  >_  5  kV  the  minimum  energy  of  the  primary  electrons 
would  exceed  the  sensitivity  band  of  the  detector  and  its  response  to  these 
electrons  would  be  zero.  Its  counting  rate  would  then  be  symmetrical  about 
a=90  since  it  would  respond  only  to  the  lower  energy  electrons  resulting 
from  interactions  in  the  atmosphere  and  local  trapping.  In  this  case  the 
potential  above  the  satellite  can  be  estimated  by  applying  the  methods 
described  above  to  the  counting  rates  of  the  CMED  detector  if  those  rates 
are  sufficiently  above  the  noise. 

3 .  RESULTS 

In  Figure  10  the  potential  difference  <j>m  -  4>s  Is  plotted  against  the 
satellite  spin  number.  The  scale  at  the  top  of  the  figure  gives  the  invariant 
latitude  and  magnetic  local  time.  Near  the  top  of  the  figure  the  shading 
indicates  the  ranges  over  which  the  field-aligned  ions  (FAI)  are  observed, 
and  the  asterisks  denote  the  more  isolated  occurrences  of  the  conical  ion 
distributions  (CID) .  The  potential  differences  are  quite  variable  and  at 
times  very  high,  exceeding  8  kV  on  one  spin  period.  Note  also  that  the 
field-aligned  ions  are  observed  whenever  e(£  -  <£s)  approximately  exceeds  the 

eneigv  threshold  of  the  ion  spectrometers. 

The  regions  in  which  the  higl  electron  fluxes  (spikes)  ,  both  near 
a=0  and  a=L80°,  are  observed  are  also  shown  in  the  figure.  Note  that  these 
spikes  generally  occur  when  <J»  -  Hss  is  very  small.  The  source  of  the  spikes 
has  not  yet  been  determined,  but  we  believe  they  may  be  due  to  torroidal, 

MHD  turbulence,  which  accelerates  electrons  in  both  directions  along  the 
magnetic  field. 

The  potential  difference  above  the  satellite  is  shown  in  Figure  11  as 
a  function  of  the  satellite  spin  number.  These  values  of  ^ s  vary  in  approxi¬ 
mately  the  same  manner  as  those  in  Figure  12.  An  important  difference,  how¬ 
ever,  is  that  the  potential  <J>S  is  always  high  just  preceding  the  occurrences 
of  the  electron  spikes.  The  magnetic  field  ratios  B^/Bs  were  found  to  be 
less  than  0.78,  with  most  of  the  separately-determined  values  being  between 
about  0.55  and  0.70. 

The  variations  of  the  potential  differences  shown  in  Figures  10  and 
11  are  too  rapid  to  be  ascribed  to  the  spatial  (altitude  and  invariant  lati¬ 
tude)  or  magnetic  local-time  variations  of  the  satellite.  As  described  in 
References  1  and  2  the  potential  differences  were  found  to  have  an  inverted- 
V  structure,  a  few  degrees  wide  in  invariant  latitude.  Furthermore,  since 
this  current  sheet  is  associated  with  the  discrete  aurorae,  rapid  variations 
are  not  expected  within  the  current  sheet  in  the  longitudinal  direction. 

The  rapid  variations  are  therefore  more  likely  due  to  oscillatory  motions  of 
the  magnetic  shell  containing  the  f ie Id -aligned  currents.  Observations  of 
oscillations,  principally  with  periods  of  2  to  7  minutes,  of  the  magnetic- 
field  at  high  latitudes  and  in  the  motion  of  the  magnetopause  suggest  that 
the  high  latitude  field  lines  fluctuate  with  such  periods.  The  variations 
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SPIN  NO. 

Figure  10.  Potential  difference  below  satellite  as  function  of  satellite  spin  number.  The 
invariant  latitude  and  magnetic  local  Ome  are  shown  at  the  top  of  the  graph.  Observations 
of  the  field-aligned  ions  ( FA  I )  and  conical  ion  distributions  (CID)  are  indicated  near  the 
top  of  the  figure.  The  locations  of  the  electron  spikes  are  also  shown  in  the  figure. 


Figure  11.  The  potential  difference  above  the  satellite  as  a  function 
of  the  satellite  spin  number.  The  arrows  indicate  that  the  <(>-  is  between 


of  the  potentials  in  Figures  10  and  11  appear  to  be  consistent  with  fluctua¬ 
tions  of  this  type. 


The  effect  of  the  variation  of  the  invariant  latitude  of  the  satellite 
is  indicated  by  the  general  lowering  of  the  activity  at  invariant  latitudes 
higher  than  about  76.2  .  The  data  described  in  References  1  and  2  display  an 
electron  spike  at  the  high  latitude  boundary  of  the  potential  structure. 

The  regions  of  the  spikes  shown  in  Figures  10  and  11,  especially  the  locations 
of  the  large  spikes  appear  to  denote  the  high  latitude  region  of  the  potential 
structure . 

Our  interpretation  of  the  results  is,  therefore,  as  follows.  On  spin 
zero  the  satellite  entered  the  region  in  which  the  narrow  potential  structure, 
due  to  the  current-driven  instability,  was  fluctuating  in  the  latitudinal 
direction.  On  spin  numbers  43-74  the  satellite  had  generally  cleared  the 
region,  on  the  high-latitude  side,  except  for  the  surge  of  the  current  sheet 
that  occurred  in  the  vicinity  of  spin  number  38.  The  satellite  then  again 
traversed  that  region  as  it  moved  toward  lower  latitudes  and  earlier  magnetic 
local  times.  The  potentials  were  generally  higher  during  the  local  evening 
than  during  the  local  afternoon. 
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SECTION  IV 

ANALYSIS  OF  ION  DISTRIBUTIONS 


1.  PITCH-ANGLE  DIFFUSION  IN  PRESENCE  OF  ELECTRIC  FIELDS 

Pitch-angle  diffusion  is  customarily  treated  with  a  Fokker-Planck 
diffusion  equation  for  the  phase  space  distribution  function,  f.  Since 
the  principal  observables  are  kinetic  energy,  w,  and  pitch-angle,  a,  f  is 
usually  regarded  as  a  function  of  w  and  a.  However,  the  proper  physical 
variables  are  the  constants  of  motion,  in  this  case  the  adiabatic  invariants. 
To  see  how  this  affects  the  form  of  the  Fokker-Planck  equation,  consider 
pitch-angle  diffusion  in  a  divergent  magnetic  field  with  energy  conservation 
(i.e.  trapped  electrons  in  a  VLF-ELF  wave  field;  see  Kennel  and  Petschek, 

Ref.  21).  The  local  interactions  are  described  by  the  homogeneous  field 
equation, 
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where  Df/Dt  is  the  total  derivative, 
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Equat ion  (1 3) cannot  represent  the  variation  of  the  distribution  along  the 
magnetic  field  direction  because  a  is  not  an  adiabatic  invariant  variable. 
The  pitch-angle  at  the  equator,  an,  is  an  adiabatic  invariant  variable, 
so  the  non-local  analog  of  Equation  (13)  is. 
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and  IL,  ~  a  bounce  averaged  diffusion  coef- 
a0a0  ^ 


In  the  presence  of  two  non-stochast ic  effects,  a  divergent  magnetic 
field  and  a  parallel  electric  field,  we  must  find  new  constants  of  the  motion. 
Conservation  of  the  first  adiabatic  invariant  (magnetic  moment)  gives  a 
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where  and  wu  represent  the  (non-relativist ic)  components  of  w.  Motion 
along  the  field  direction  must  preserve  the  total  energy,  kinetic  plus  poten¬ 
tial.  If  v  is  the  total  potential  energy  (charge  e  times  electric  potential 
M  we  have  the  definition  of  another  constant  C^: 
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C2  -  CjB 


Wjj  =  w(l  -  sin  a) 


the  zero  point  for  *0  is  arbitrary,  so  we  can  choose  the  equator  to  be  the 
point  where  tf;  =  0.  This  alows  us  to  construct,  from  the  above  equations, 
two  new  constants  of  motion. 


u  =  i  -  ^  (l£ 

K  -  w  +  i'  (1$ 

where  BQ  is  the  magnetic  field  at  the  equator  (or  ^  =  0  point).  Clearly  U 
is  related  to  a0,  for  it  becomes  equal  to  cos2  aQ  at  the  equator;  i.e. 
the  second  term  on  the  right  hand  side  of  (18)  becomes  w  / w,  which  is  equal 
to  s in^  cxq  . 

The  derivative  operator  of  Equation  13  is 
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Substituting  this  in  the  diffusion  equation  gives 
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where  f  must  be  a  function  of  w  -  K  and  U.  Note,  however  that  (21)  is  still 
a  local  diffusion  equation. 

To  obtain  the  bounce-averaged  equation  integrate  over  the  particle’s 
trajectory.  The  element  of  length  is  ds/cos  ,  where  ds  is  the  increment 
along  the  magnetic  field;  the  integral  is 
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where  and  S2  are  the  locations  of  the  conjugate  mirror  points,  and  T  is 
the  quarter- bounce  integral. 
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R0  is  the  distance  from  the  center  of  the  earth  to  the  equatorial  crossing  of 
the  magnetic  field  line. 

The  development  of  Equation  (22)  is  the  same  as  for  Equation  (15),  as 
can  be  seen  immediately  by  letting  ip  ~  0  and  U  -  cos^  aQ.  The  left  side  of 
the  equation  is  not  so  simple,  however,  because  the  electric  field  can  raise 
or  lower  the  mirror  point.  If  the  mirror  points  are  both  above  the  atmosphere, 
the  effects  of  absorption  can  be  neglected  and  the  ^  .  Vf  term  is  near  zero. 

If  either  mirror  point  dips  into  the  atmosphere,  the  absorption  contributes 
a  term  of  order 
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The  diffusion  equation  above  is  appropriate  to  interactions  with 
electromagnetic  (cyclotron-mode)  waves  that  approximately  preserve  the 
particle  energy.  Other  waves,  such  as  transverse  electrostatic  waves,  prim¬ 
arily  alter  the  perpendicular  component  of  momentum.  The  local  equation  is 
then. 
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In  this  case  there  is  no  simple  decomposition  into  an  adiabatic  invariant 
equation  with  a  single  degree  ot  freedom  (as  in  (15)  or  (22)).  This  can  be 
seen  most  readily  by  finding  dll/dK  and  noting  that  it  is  not  a  single-valued 
function  of  U  and  K.  The  rightside  must  comprise  several  terms,  including 
cross  terms  with  two  different  derivatives.  It  is  highly  unlikely  that  analy¬ 
tic  solutions  exist;  even  numerical  solutions  are  difficult  in  an  inhomo¬ 
geneous  magnetic  field. 


A  useful  set  of  variables  is  K  and 
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This  is  simply  the  first  adiabatic  invariant,  in  energy  units.  The  differen¬ 
tial  operator  is 


(29) 


A3 


y,v 


The  diffusion  equation  becomes 


Df  _  (ji  +  _3_\  f  (; 

Dt  \  B  3M  +  3K/  |DJ-L\ 


_o  _3_  _3_ ' 

B  5M  +  3K  J 


K  8D^  .  !i±\  Km  .  at) 

\  B  3M  3K  /  \  B  3M  3K/ 


/B  .2_  B  .2,  .2. 

D  f  — - +  2  —  — f 

X1\b2  3M2  B  9M8K  3K2 


The  cross  terms  cannot,  of  course,  be  eliminated,  so  the  equation  can  have 
all  combinations  of  first  and  second  derivatives  on  the  right. 

2.  DIFFUSION  COEFFICIENTS  FOR  INTERACTIONS  WITH  ELECTROSTATIC  WAVES 

The  Fokker-Planck  equation  for  interactions  that  alter  the  perpendicular 
and  parallel  parts  of  a  particlefs  energy  is  (Chandresekhar ,  Ref.  23), 
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The  notation  <AX>  for  the  Fokker-Planck  coefficients  denotes  a  time  average 
of  the  changes  in  X,  per  unit  time.  If  the  waves  affect  only  one  component, 

Aw  ,  the  w  terms  disappear.  An  equilibrium  solution  of  the  remaining  equa¬ 
tion  is  3f/3t  =  0,  f  =  constant.  This  gives  a  relation  between  the  Fokker- 
Planck  coefficents, 
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A  solution  of  (12)  is 
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The  Fokker-Planck  equation  can  then  be  written  as  a  diffusion  equation 
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Interactions  of  charged  particles  with  transverse  electrostatic  waves 
may  be  assumed  to  affect  only  Wj_#  A  particle  of  charge  e  in  resonance  with 
a  transverse  wave  field,  of  amplitude  Ew  experiences  a  change  in  energy. 
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If  the  wave  and  particle  are  out  of  resonance  by  an  amount  6w,  they  will  re¬ 
main  in  phase  for  a  time  of  order 
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Generally,  except  in  a  highly  inhomogeneous  field,  the  interaction  time  is  set 
by  relation  (37) .  The  diffusion  coefficient  follows  upon  taking  a  time  average 
of  (Aw^)2,  thus, 
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The  last  term  on  the  right  has  the  form  of  a  power  spectrum.  The  average  wave 
power  in  a  band  of  width  6v  =  6 os/2  is  the  power  spectral  density, 
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The  diffusion  coefficient  therefore  becomes, 
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(40) 


The  power  spectral  density  of  the  broad-band  electrostatic  turbulence,  in 
the  region  of  the  peak  intensity  (10  -  50  Hz),  measured  by  Gurnett  and  Frank 
(Ref.  9)  is  about  5  x  10“^  (V/m)^/Hz.  Substituting  this  value  in  (40)  gives, 
for  protons, 

D  =  0.06  w  (keV)^/sec  (41) 

Here,  the  energy  is  in  keV.  A  more  detailed  analysis  by  Sturrock  (Ref.  24) 
led  to  a  similar  form,  with  proportional  to  w^. 
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3. 


SOLVING  THE  DIFFUSION  EQUATION 


The  local  diffusion  equation,  (35),  has  the  form. 


11  =  r  L  1L\ 

3t  3wv  \  3-  3W  / 


The  initial  condition  is  that  the  distribution  is  a  delta  function  at  =  0. 
The  boundary  condition  at  =  0  is  w^  8f/3w_^  =0.  A  solution  that  satisfies 
these  conditions  is 

f  ~  ni  exp  <-wx/rt)  (43) 

This  solution  is  useful  for  verifying  the  computational  results  at  late  time, 
but,  unfortunately,  the  adiabatic  invariant  diffusion  equation  has  no  such 
simple  solutions. 

An  alternative  that  avoids  the  difficulty  of  solving  the  diffusion  equa¬ 
tion  is  to  return  to  the  fundamental  processes  and  trace  the  probable  trajec¬ 
tories  of  individual  particles.  In  our  case  a  Monte  Carlo  technique  is  well 
suited  to  individual  trajectories  if  the  concept  of  mean  deflection  is  invoked. 
We  postulate  that  the  deflection,  6w^  ,  obeys  a  probability  distribution  (nor¬ 
malized  to  1)  :  ~ 


P(6w,)  = 


exp  [ - ( <$w,  )^/2 a^] 


This  distribution  is  not  based  on  physical  principles,  but  rather  has  the 
expected  form  and  a  well-defined  mean,  cj2;  there  are  few  processes  in  nature 
that  strictly  obey  a  normal  (Gaussian)  distribution.  The  critical  restric¬ 
tions  on  the  use  of  the  probability  distribution,  (44) ,  are  that  the  mean 
deflection  be  small  compared  with  the  scale,  f/(df/dw(),  of  the  distribution 
and  small  compared  with  Wj_.  These  conditions  can  be  "satisfied  by  choosing 
the  time  interval  to  be  sufficiently  small. 

The  Fokker-Planck  coefficients  are  related  to  P.  The  diffusion  coef¬ 
ficient  is 


,  (6w,r 
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2At  /  P(5wl)(5wi.)  d{wL  =  2At  =  W-Lr  (A5) 


The  probability  distribution  is  now 


P  ( 6wx) 
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exp  ( -(6w. ) / (4w.  FAt)) 


A  random  variable,  Q,  uniformly  distributed  between  0  and  1,  was  determined  from 
the  cumulative  probability  distribution. 


Q(6w. )  =  /  P(x)  dx 


At  the  end  of  each  predetermined  time  step,  At,  Q  was  set  equal  to  a  random 
number  and  (47)  was  solved  for  5w^_.  Negative  as  well  as  positive  values  of 
6wj_  were  allowed.  Negative  values  of  w^  of  course  are  not  physically  possible. 

Two  methods  were  used  to  deal  with  this  computational  possibility.  In  one,  the  ions 

were  assumed  to  be  '’absorbed”  at  w^  <  0,  i.e.  the  case  history  was  eliminated 

if  w  became  negative.  In  the  other,  reflection  about  Wj_  =  0  was  assumed; 

the  sign  of  w^  was  reversed  if  it  became  negative.  The  ion  distributions 

computed  on  the  basis  of  these  widely  different  assumptions  were  essentially 

the  same,  undoubtedly  because  the  diffusion  coefficient  favored  the  higher 

values  of  w^. 

The  effect  of  the  divergence  of  the  magnetic  field  on  the  ion  trajectories 
was  taken  into  account  by  assuming  that,  between  resonances,  the  adiabatic 
motions  of  the  ions  held.  The  following  equations  were  used  to  compute  the 
ion  trajectories: 
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The  program  began  by  computing  the  trajectories  of  ions,  with  specified 
initial  values  of  wjj  and  Wj_,  from  an  altitude  Hso,  at  the  base  of  the  source 
distribution,  toward  higher  altitudes  along  the  magnetic  field.  At  certain 
"read-out"  altitudes,  Hj^,  above  Hgo,  the  characteristics  of  the  ions  were 
tabulated:  the  ions  were  accumulated  in  energy  cells  of  widths  (Aw) ^  and  in 

pitch-angle  cells  of  solid-angle  widths  (AP)^.  Also,  the  ions  in  each  of  the 
energy  cells  were  further  subdivided  into  pitch-angle  cells  of  solid-angle 
widths  (A £0^.  After  the  trajectories  of  NQ  ions,  starting  from  Hsq  were 
computed  and  tabulated  in  this  manner,  the  trajectories  of  N0  exp [ -AHs/Hscaj_e ] 
ions,  starting  from  the  altitude  Hso  +  AHS,  were  computed,  continuing  the 
tabulation  described  above.  Here,  AHS  is  the  increment  in  the  source  altitude, 


.t  *»,*** 


and  Hsca^e  is  the  scale  height  of  the  source.  The  computation  continued,  with 
ion  trajectories  starting  at  ever  higher  altitudes,  until  the  source  altitude 
exceeded  the  maximum  readout  altitude. 


The  incremental  directional  flux  at  a  readout  altitude  due  to  the  ions 

AN  (s  )  that  originated  at  the  altitude  s  within  the  volume  element  A  AH 
,  ik  n  n  ns 

is 


ANik(Sn} 

Aj(wi?  a^,  sn)  Aw.Aft,  Acosa,  At 


(53) 


where  A  is  the  cross  sectional  area  of  the  magnetic  tube  at  the  readout 
altitude.  But,  from  the  continuity  equation, 
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where  vOJl  is  the  initial  parallel  velocity  component  of  the  ions.  By  putting 
(34)  into  (53)  and  summing  over  the  source  altitudes  s  ,  the  flux  is  given  by 
the  equation,  n 


V  V  AN.,  (s  ) 
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i’  k  Aw , Aft.  Acosa.  AH 
i  k  k  s 


(55) 


If  AHS  would  have  been  variable  along  the  tube,  the  summation  in  (55)  would 

have  been  over  the  ratio  AN.,  (s  )/AH 

ik  n  sn 

The  following  information  was  printed  at  each  readout  altitude: 


(i)  the  number  of  ions  in  the  energy  cells,  N.,  and  in  the  solid- 

angle  cells,  N^,  1 

(ii)  the  number  of  ions,  N.,  in  the  energy  cells  divided  by  the 
widths  Aw^  (this  ratio  is  proportional  to  the  energy  distri¬ 
bution  of  the  flux,  i.e. ,  to  the  summation  of  (55)  over  the 
k  variables) 
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(iii)  the  number  of  ions  Nik  with  energies  in  the^ith  cell  and  pitch- 
angles  in  the  kth  cell  divided  by  AWj^A^^cosa^,  where  ak  is  the 
mean  value  of  a  in  the  kth  cell  (this  ratio  is  proportional  to 

j  as  given  by  Eq .  55),  and, 

(iv)  the  number  of  ions  Nk  in  the  pitch- angle  cells,  for  all  energies, 
divided  by  Aft,  cosa,  (this  ratio  is  proportional  to  the  pitch- 
angle  distribution  of  the  ions,  i.e.,  to  the  summation  of  (55) 
over  the  ith  variables. 


4.  RESULTS  OF  MONTE  CARLO  CALCULATION 

The  Monte  Carlo  program  was  first  run  with  E  =0  and  VB  =  0  to 
determine  whether  the  numerical  solution  of  the  diffusion  equation  (42)  would 
compare  favorably  with  the  analytical  solution  (43).  The  following  parameters 
were  used:  heating  rate,  V  =.0283  keV/sec;  diffusion  time,  t  =  27.95  sec; 

At  =  .2795  sec  (100  collisions);  initial  energy,  =  .5  keV ;  and  number  of 
case  histories,  N  =  500.  The  resulting  energy  distribution,  f,  of  the  ions 
is  shown  in  Figure  12.  The  straight  line  drawn  in  this  figure  has  the  slope 
Ft  *  .79  keV,  which  is  the  form  of  the  analytical  solution.  Note  that  the 
computed  values  (the  open  circles)  follow  the  analytical  solution  quite  well. 
Actually,  the  computed  values  are  still  somewhat  elevated  in  the  vicinity  of 
.5  keV,  the  initial  value  of  w±.  Hence,  a  shorter  value  of  At  (more  collisions 
within  the  time  t)  would  be  desirable. 

Note  from  (41)  and  (45)  that  the  heating  rate,  r ,  used  in  this  run  is 
comparable  to  the  value  .06  keV/sec  implied  by  the  measurements. 

A  series  of  runs  was  then  conducted,  with  Ejj  >  0  and  VB  appropriate  for 
the  earth's  field,  to  determine  the  effects  of  the  model  parameters  on  the 
energy  and  pitch-angle  distributions  of  the  ions.  Combinations  of  the  para¬ 
meters  were  sought  that  would  give  distributions  that  are  representative  of  the 
measured  distributions.  In  references  1  and  2,  three  measurements  of  the  energy 
distributions  of  the  ion  flux  made  during  the  S3-3  satellite  traversal  of  the 
inverted-V  potential  structure  are  described.  These  energy  spectra  are 
reproduced  in  Figure  13.  The  spectra  in  the  panels  from  left  to  right  will  be 
referred  to  by  the  numbers  1,  2,  and  3,  respectively.  At  the  times  of  the 
measurements  of  the  spectra  1,  2,  and  3,  the  potential  differences  below  the 
satellite,  inferred  from  the  electron  distributions,  were  about  .8,  2,  and 
2.5  keV,  respectively.  Using  these  potentials  and  typical  observed  wave-power 


Figure  13.  Energy  spectra  of  ions  measured  by  the  Lockheed  and  Aerospace 
groups  on  the  three  consecutive  spins  of  the  S3-3  satellite  through  the 
inverted-V  structure  discussed  in  References  1  and  2. 
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spectra  as  constraints  on  the  model  parameters,  the  combinations  of  the  para¬ 
meters  that  resulted  in  spectra  similar  to  spectra  1  and  2  are  listed  in 
Table  2. 


Table  2.  Parameters  in  model  that  give 
computed  ion  spectra  similar  to  observed  spectra  1  and  2. 


Spectrum 

H 

so 

H 

scale 

hr 

En 

r 

At 

Number 

(km) 

(km) 

(km) 

(mV/m) 

(keV/sec ) 

(sec) 

1 

4000 

600 

8000 

0.2 

.0306 

.04083 

2 

4000 

600 

8000 

0.5 

.1061 

.06455 

The  corresponding  pitch-angle  distributions  were  consistent  with  the  measure¬ 
ments,  but  these  measurements  were  not  sufficiently  precise  to  provide  an 
adequate  test  of  the  computational  results. 

Some  of  the  computed  ion  distributions  are  shown  in  Figures  14-29. 

All  these  distributions  were  computed  using  the  following  parameters: 

Hso  =  4000  km,  AHS  =  400  km,  Hsca^e  -  600  km,  w^q  =  .01  keV,  and  wxq  =  .3  keV. 
The  number  of  case  histories  from  the  base  altitude,  NQ ,  was  2000.  The  values 
of  E  || ,  F,  At,  Hr,  and  A(f>  =  E|,(Hr  =  Hso)  are  given  in  the  figure  captions.  The 
histograms  shown  in  Figures  14  and  15  are  the  energy  spectrum  of  the  flux 
(print-out  item  (ii)  described  above)  and  the  corresponding  pitch-angle 
distribution  of  the  ions  (print-out  item  (iv))  at  Hr  =  6000  km,  for  E ^  =  .2  mV/m, 
A(j>  =  .4  kV,  T  =  .0433  keV/sec,  and  At  =.02887  sec.  The  ordinates  are  related 
to  the  counts  in  the  cells  as  described  in  items  (ii)  and  (iv)  above.  Note 
that  both  distributions  are  quite  broad,  owing  to  the  high  energy  that  the 
ions  acquire  from  the  transverse  heating  mechanism  relative  to  the  energy  that 
the  ions  acquire  by  falling  through  the  potential  drop  along  the  magnetic 
field.  The  output  information  described  in  item  (iii)  above  reveals  that  the 
pitch-angle  distributions  of  the  ions  are  strongly  energy  dependent;  toward 
increasing  ion  energies,  the  pitch-angle  distributions  become  narrower  and 
peaked  at  larger  pitch  angles.  Hence,  conical  pitch  angle  distributions  may 
be  observed  not  only  when  the  ions  are  rapidly  accelerated  within  a  short 
distance  along  the  magnetic  field  by  coherent  EIC  waves,  but  also  when  the 
ions  are  accelerated  more  slowly  by  the  waves  over  a  long  distance  if  the 
energy  threshold  of  the  detector  is  sufficiently  high;  in  this  case,  about 
1.3  keV. 

The  distributions  at  the  higher  altitude,  Hr  =  8000  km  computed  with  the 
parameters  listed  for  spectrum  number  1  in  Table  2,  are  shown  in  Figure  16  and 
17.  The  measured  spectrum  is  also  drawn  in  Figure  16.  Note  that  the  computed 
distribution  is  in  fair  agreement  with  the  measured  value.  The  number  density 
of  the  ions  in  Spectrum  1  is  about  0.2/cm-^  (Ref.  1).  Application  of  the 
continuity  equation  reveals  that  if  all  these  ions  originated  4000  km  below 
the  satellite  with  initial  energies  of  .01  keV,  the  number  density  of  the  ions 
at  the  source  would  be  about  5/cm^. 
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Figure  18.  Histogram  of  directional  flux  versus  energy 
E ||  =  .2  mV/n,  Y  =  .08660  keV/sec,  At  =  .02887  sec, 

=  60r.0  km,  and  =  .4  kV. 


Figure  28.  Histogram  of  directional  flux  versus  energy 
K  =  .5  nV/m,  p  =  .1061  koV/sec,  At  =  .06455  see, 

0R  -  8000  km,  and  =  2  kV.  The  curve  is  the  measured 
flux  (spectrum  number  2)  fitted  to  the  histogram. 
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A  comparison  of  the  distributions  in  Figures  14  and  16  indicate  that 
Spectrum  1  could  also  have  been  matched  if  the  base  of  the  source  had  been 
chosen  to  be  2000  km  (i.e.,  HR  -  Hso)  below  the  satellite  and  E  had  been 
increased  to  0.4  mV/m. 

The  succeeding  figures,  18-25,  principally  illustrate  the  effect  of 
an  increasing  heating  rate  on  the  distributions  at  the  two  altitudes,  6000  km 
and  8000,  above  the  source  base,  4000  km.  The  effect  of  an  increased  electric 
field  is  illustrated  in  Figures  26-29.  Spectrum  2  is  compared  with  the  ion 
distribution  in  Figure  28.  Here,  again,  the  agreement  is  seen  to  be  fairly 
good.  It  could  have  been  improved  by  decreasing  somewhat  and  increasing  T. 
Spectrum  2  could  also  have  been  fit  by  choosing  parameters  similar  to  those 
listed  in  the  caption  of  Figure  26  (a  source  base  2000  km  below  the  satellite) 
but  doubling  the  electric  field,  to  1  raV/m. 

The  distribution  in  Figure  26  matches  very  well  the  measured  spectrum, 

3.  However,  in  the  computation  of  that  distribution,  A<j>  was  0.8  kV.  The 
potential  difference  below  the  satellite  at  the  time  of  the  Spectrum  3  measure¬ 
ment  was  about  2.5  kV.  In  the  context  of  the  model  discussed  above,  a  peak 
in  the  spectrum  below  the  energy  eA can  only  be  explained  if  the  source  inten¬ 
sity  increases  at  altitudes  above  the  potential  4>m.  The  potential  difference 
above  the  satellite  was  very  small  when  spectra  1  and  2  were  measured;  it  then 
increased  sharply  on  the  following  satellite  spin,  when  spectrum  3  was  measured, 
implying  an  intensification  of  the  wave  field  above  the  satellite.  If  the 
wave  field  also  increased  below  the  satellite,  but  above  the  point  where  <J>  =  4>m, 
an  increase  in  the  ion  source  intensity  above  the  base  could  be  explained. 


SECTION  V 


CONCLUSIONS 


The  data  obtained  on  17  January  1977,  while  the  S3-3  satellite  was  along 
a  long  trajectory,  from  21. 3h  MLT  to  16. Oh  MLT,  through  the  nominal  region  of 
the  upward  magnetic-field-aligned  current  sheet,  yielded  important  new  informa¬ 
tion  on  the  properties  of  the  potential  structure.  The  results  described  in 
Section  III. 3  indicate  that  high  potential  differences,  to  values  exceeding 
10  kV,  along  the  magnetic  field  were  present  over  the  entire  local-time  interval 
of  the  satellite  pass  through  the  region.  The  results  also  indicate  that  ener¬ 
getic  H+  and  0+  ions  with  pitch-angle  distributions  generally  peaked  along  the 
magnetic  field,  but  sometimes  peaked  at  angles  inclined  to  the  field  (conics), 
were  present  whenever  the  potential  differences  occurred. 

Ghielmetti  et  al.  (Ref.  4),  by  examining  the  data  on  the  upward-flowing 
ions  obtained  with  the  S3-3  satellite  on  traversals  of  the  current-sheet  region 
at  different  local  times  but  within  short  local-time  intervals,  found  that  the 
probability  of  observing  upward-flowing  ions  over  the  auroral  oval  is  high 
(>60%  above  an  altitude  of  4000  km).  The  probability  of  such  observations  is 
the  highest  during  the  local  evening.  At  other  local  times,  the  probability 
is  highly  correlated  with  magnetic  activity. 

The  results  of  this  statistical  study,  therefore,  together  with  the 
observations  in  the’ current  sheet  region  over  the  wide  local-time  interval, 
indicate  that  the  current  sheet  on  the  dusk  half  of  the  auroral  oval  is 
nearly  always  unstable  over  a  wide  range  of  longitude.  This  finding  is  con¬ 
sistent  with  the  model  of  the  antisolar  convection  of  the  polar-cap  field 
lines  over  the  polar  cap  discussed  in  Section  II  and  with  the  instability 
criteria  provided  by  Kindel  and  Kennel  (Ref.  14).  Potential  difference  along 
the  magnetic  field  evidently  appear  on  the  dawn  side  of  the  auroral  oval  at 
high  values  of  Kp  when  the  region  2  current  system  develops. 

The  latitudinal  fluctuations  of  the  potential  region  inferred  from  the 
data  (see,  e.g.  Figure  10)  are  not  surprising  in  view  of  the  behavior  of  the 
discrete  aurorae  and  the  expected  resonant  oscillations  of  high-latitude 
field  lines.  An  important  consequence  of  this  finding,  however,  is  that 
latitudinal  motions  of  the  current  sheet  induce  electric  fields  along  the 
magnetic  field.  The  nature  and  effects  of  such  electromagnetic  fields  have 
not  yet  been  assessed. 

Another  fascinating  feature  of  the  data  are  the  spikes — the  high  fluxes 
of  electrons  that  are  highly  collimated  in  both  directions  along  the  magnetic 
field.  The  source  of  these  spikes  has  not  yet  been  investigated. 


The  analysis  of  the  ion  distributions  indicate  that  the  ions  are  simul¬ 
taneously  heated  transversely  by  interacting  with  the  EIC  wave  turbulence  and 


accelerated  longitudinally  by  the  electric  field  along  the  magnetic  field. 

The  simulation  of  this  acceleration  process  in  a  Monte  Carlo  program,  using 
a  transverse  heating  rate  appropriate  for  typical  values  of  wave  power  spectra 
and  potential  differences  below  the  satellite,  was  successful  in  reproducing 
the  usual  characteristics  of  the  ions  in  the  potential  region.  The  results 
indicate  that  energy  spectra  with  peaks  at  energies  equal  to  or  less  than  the 
potential  energy  below  the  satellite  can  be  reconstructed  by  assuming  the  ion 
source  to  be  extended  along  the  magnetic  field,  with  its  base  located  at 
several  thousand  kilometers  below  the  satellite  and  its  intensity  diminishing 
exponentially  toward  higher  altitudes  with  a  scale  height  less  than  about  1000 
km.  If  the  peaks  of  the  energy  spectra  are  at  energies  less  than  the  potential 
energy  below  the  satellite,  the  source  intensity  must  increase  with  altitude. 

In  all  the  simulation  runs,  the  electric  field  was  taken  to  be  equal  to  the 
potential  difference  below  the  satellite  divided  by  the  distance  between  the 
source  base  and  the  satellite. 

The  pitch-angle  distributions  of  the  ions  resulting  from  this  accelera¬ 
tion  process  are  both  peaked  at  a  =  180  and  conically-shaped,  depending  on 
the  energy  band  of  the  ions.  At  low  energies,  the  distributions  are  peaked 
along  the  field;  toward  increasing  ion  energies,  the  distributions  become 
peaked  at  increasing  inclination  angles. 

The  distribution  of  the  source  intensity  used  in  the  simulation  program 
may  reflect  solely  the  intensity  of  the  wave  field,  since  a  source  number 
density  less  than  about  5  ions/cm^  is  needed  to  account  for  the  observed  ion 
fluxes.  Hence,  the  distribution  of  the  wave  turbulence,  as  well  as  the  scale 
of  the  potential  distribution,  are  found  to  be  consistent  with  the  anomalous- 
resistivity  theory. 


SECTION  VI 


RECOMMENDATIONS 


The  close  relationship  between  the  phenomena  occurring  in  the  polar- 
cap  regions  and  in  magnetic  tubes  containing  hot  plasma  resulting  from  a 
nuclear  explosion  has  been  established.  The  instability  criteria  predicted 
by  Kindel  and  Kennell  (Ref.  14),  and  verified  in  the  auroral  current  sheet, 
indicate  that  the  currents  at  the  boundaries  of  tubes  containing  hot  plasma 
due  to  megaton-type  explosions  will  become  unstable  and  thereby  impair  the 
coupling  of  the  ionosphere  and  magnetosphere.  In  order  to  predict  the  ensuing 
effects,  more  detailed  information  on  the  natural  phenomenon  should  be  obtained. 

Several  recommendations  stem  from  this  work.  Recently,  the  S3-3  satel¬ 
lite  has  made  several  passes  through  the  region  of  the  current  sheet,  in  the 
southern  hemisphere,  similar  to  the  one  on  17  January  1977.  Since  the  analysis 
of  the  17  January  data  has  yielded  so  much  information  on  the  properties  and 
dynamics  of  that  region,  analysis  of  further  data  of  this  type  is  recommended. 
Attempts  should  be  made  to  correlate  such  observations  with  simultaneous  measure¬ 
ments  in  or  near  the  current  sheet  of  the  wave  turbulence  and  ion  measurements 
on  the  GEOS  and  ISEE  satellites. 

fhe  high  fluxes  of  counterstreaming  electrons  (the  electron  spikes) 
continue  to  pose  a  perplexing  and  fascinating  problem.  Further  data  on  this 
phenomenon  should  be  analyzed,  and  an  intensive  theoretical  investigation 
should  be  undertaken.  Such  a  study  may  very  likely  uncover  a  dynamical  plasma 
process  which  may  have  an  important  impact  on  some  of  the  defense  problems 
discussed  in  Section  I. 

The  electric  field  induced  by  the  motion  of  the  current  sheet  above  the 
auroral  oval,  as  well  as  the  effects  of  this  electric  field,  should  also  be 
investigated. 

Information  on  the  location  and  extent  of  the  region  where  the  ions  are 
accelerated  perpendicularly  to  the  field  can  be  inferred  fairly  directly  from  the 
energy  and  conical  pitch-angle  distributions  of  the  ions.  Wide  pitch-angle 
distributions  at  high  ion  energies  indicate  a  large  extent  of  the  heating 
region.  The  Monte  Carlo  program  could  be  used  for  the  analysis  of  such  cases. 
Maxwellian  energy  distributions  and  narrow  pitch-angle  distributions  with  peaks 
at  angles  that  increase  with  increasing  energy  indicate  a  short  extent  of  the 
heating  region.  Such  cases  could  be  investigated  analytically  through  an  appli¬ 
cation  of  Liouville's  theorem.  Figure  30  depicts  the  form  of  the  flux,  as  a 
function  of  pitch  angle,  expected  for  ions  that  have  been  heated  transversely 
to  a  temperature  Tt  -  0.8  keV  below  the  satellite,  where  the  magnetic  field 
intensity  was  5  times  the  field  intensity  at  the  satellite  and  where  the 
potential  was  1  kV  higher.  The  arrows  denote  the  pitch-angles  of  the  ions 
which  have  the  energies  shown  at  the  arrows.  The  ion  temperature  and  number 
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igure  30.  Directional  flux  at  Bs ,  versus  pitch  angle  for  ions  that  have  been 

ransversely  heated  at  B  =  5B  ,  y  -  =  1  kV  to  a  temperature  of  0.8  keV.  The 

oints  noted  on  the  curve  are  the  pitch  angles  of  the  ions  which  have  the  energies 
hown  at  the  arrows. 


density,  and  the  location  and  potential  of  the  source  point,  can  be  inferred 
from  such  distributions.  Finally,  narrow  energy  distributions,  in  addition 


to  narrow  pitch-angle  distributions,  imply  acceleration  by  coherent  EIC  waves. 
Such  distributions  can  also  be  analyzed  by  using  Liouville's  theorem.  These 
analyses  require  highly  resolved  measurements  of  both  the  energy  and  pitch- 
angle  distributions  of  the  ions,  in  addition  to  the  electron  measurements 
which  are  required  for  the  determination  of  the  potential  difference.  Approp¬ 
riate  measurements  are  available  from  the  combined  experiments  of  Lockheed 
and  Aerospace  on  the  S3-3  satellite. 

In  addition  to  the  investigations  of  the  natural  phenomena  recommended 
above,  the  possible  effects  of  anomalous  resistivity  in  magnetic  tubes  con¬ 
taining  hot  plasma  produced  by  nuclear  bursts  should  be  estimated.  For  these 
investigations,  the  potential  difference  along  the  magnetic  field  should  be 
assumed  to  be  a  sizable  fraction  of  the  potential  difference  across  the  tube, 
say  l/4th,  as  in  the  case  of  the  polar-cap  region. 

The  effects  of  an  electric  field  along  the  magnetic  field  and  transverse 
ElC-wave  heating  on  the  distributions  of  the  debris  in  the  magnetosphere  should 
also  be  estimated.  If,  owing  to  the  transverse  heating,  an  appreciable  fraction 
of  the  debris  will  become  trapped  by  the  magnetic  field,  the  intensity  of  the 
trapped  electrons  will  be  appreciably  enhanced.  Moreover,  the  build-up  of 
trapped  debris  in  the  radiation  belt  in  the  event  of  a  nuclear  exchange  would 
present  a  further  hazard  to  satellite  systems. 
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